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Experiments  in  a  rapid  compression  machine  were  used  to  examine  the  influences  of  variations  in  fuel 
molecular  structure  on  the  autoignition  of  isomers  of  pentane.  Autoignition  of  stoichiometric  mixtures  of 
the  three  isomers  of  pentane  were  studied  at  compressed  gas  initial  temperatures  between  640  K  and  900 
K  and  at  precompression  pressures  of  300  and  400  torr.  Numerical  simulations  of  the  same  experiments 
were  carried  out  using  a  detailed  chemical  kinetic  reaction  mechanism.  The  results  are  interpreted  in  terms 
of  a  low-temperature  oxidation  mechanism  involving  addition  of  molecular  oxygen  to  alkyl  and  hydroper- 
oxyalkyl  radicals.  Results  indicate  that  in  most  cases,  the  reactive  gases  experience  a  two-stage  autoignition. 
The  first  stage  follows  a  low-temperature  alkylperoxy  radical  isomerization  pathway  that  is  effectively 
quenched  when  the  temperature  reaches  a  level  where  dissociation  reactions  of  alkylperoxy  and  hydro- 
peroxyalkylperoxy  radicals  are  more  rapid  than  the  reverse  addition  steps.  The  second  stage  is  controlled 
by  the  onset  of  dissociation  of  hydrogen  peroxide.  At  the  highest  compression  temperatures  achieved,  little 
or  no  first-stage  ignition  is  observed.  Particular  attention  is  given  to  the  influence  of  heat  transfer  and  the 
importance  of  regions  of  variable  temperature  within  the  compressed  gas  volume.  Implications  of  this 
work  on  practical  ignition  problems  are  discussed. 


Introduction 

Hydrocarbon  ignition  is  important  in  many  prac¬ 
tical  combustion  systems,  including  internal  com¬ 
bustion  engines,  detonations,  pulse  combustors,  and 
flame  initiation.  The  rapid  compression  machine 
(RCM)  is  used  frequently  to  study  kinetics  of  hydro¬ 
carbon  autoignition  [1-8],  since  reactive  gas  tem¬ 
peratures  are  similar  to  those  in  automotive  engines 
during  diesel  ignition  and  end-gas  autoignition  in 
spark-ignition  engines. 

The  RCM  provides  a  rich  environment  for  study 
of  hydrocarbon  oxidation,  including  degenerate 
chain  branching,  alkylperoxy  radical  isomerization, 
and  thermal  feedback  [9].  Hydrocarbon  oxidation 
studies  in  the  RCM  have  been  summarized  [10,11], 
and  many  classes  of  fuels  have  been  studied.  De¬ 
tailed  kinetic  modeling  is  another  tool  to  study  au¬ 
toignition  in  the  RCM  [4,12],  The  present  work  was 
intended  to  determine  experimentally  the  influence 
of  variations  in  fuel  molecular  structure  on  autoig¬ 
nition  and  to  use  a  kinetic  model  to  understand  those 
variations.  This  study  is  unique  because  although 
other  studies  have  addressed  effects  of  variations  in 
pressure  and  equivalence  ratio  [12],  this  work  ad¬ 
dresses  effects  of  variations  in  fuel  molecular  struc¬ 
ture  for  all  isomers  of  a  single  fuel,  pentane,  in  an 
RCM.  The  three  isomers  of  pentane  possess  many 


of  the  structural  elements  that  determine  such  au¬ 
toignition  characteristics  as  octane  number  and  vari¬ 
ability  in  cool  flame  production. 

Experimental  Results 

The  Lille  experimental  facility  rapidly  compresses 
reactive  mixtures  of  hydrocarbon/oxygen/inert  by  a 
piston  stroke.  The  overall  compression  ratio  was 
9.28,  and  compression  time  in  this  study  was  60  ms. 
The  piston  stops  at  the  end  of  its  travel  (top  dead 
center,  TDC)  where  the  reacting  mixture  remains  at 
constant  volume.  Pressure-time  data  during  the 
compression  and  constant  volume  periods  were 
measured  by  a  Kistler  transducer  and  recorded  at  a 
sampling  rate  of  25  kHz.  Light  emission  associated 
with  the  cool  flame  was  measured  by  a  photomulti¬ 
plier  fitted  with  a  blue-colored  filter.  Full  details  of 
the  apparatus  and  experimental  procedures  have 
been  described  [13,14], 

Pressure  and  temperature  at  TDC  must  be  known 
as  precisely  as  possible.  Pressure  at  TDC  is  mea¬ 
sured  experimentally,  but  temperature  at  TDC  is  not 
readily  accessible  by  measurement  and  must  be  cal¬ 
culated  from  the  pressure  according  to  a  suitable 
model.  The  assumption  of  adiabatic  compression  is 
unrealistic  in  the  present  case  because  at  least  two 
factors  can  alter  the  temperatures  during  compres¬ 
sion,  including  heat  losses  to  the  chamber  walls  and 
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heat  release  from  chemical  reactions.  Experimental 
pressures  at  TDC  are  usually  lower  than  theoretical 
pressures  calculated  with  the  adiabatic  law,  suggest¬ 
ing  that  heat  losses  during  compression  are  greater 
than  reactive  heat  release. 

Therefore,  temperatures  at  the  end  of  compres¬ 
sion  were  calculated  from  the  pressure  at  TDC,  the 
initial  pressure  and  temperature,  the  compression 
ratio,  and  the  perfect  gas  law.  Temperatures  com¬ 
puted  in  this  manner  represent  spatial  averages  be¬ 
cause,  due  primarily  to  heat  losses  to  the  combustion 
chamber  walls,  different  regions  in  the  compressed 
gas  have  local  temperatures  that  are  different  from 
the  averaged  temperature.  Autoignition  delay  times 
computed  from  this  average  temperature  can  there¬ 
fore  produce  results  either  longer  or  even  shorter 
than  the  experimentally  determined  values  [12], 

The  reaction  mixture  at  TDC  is  not  homogeneous 
but  includes  pockets  of  hotter  gas  that  do  not  ex¬ 
perience  heat  loss,  as  observed  by  laser-induced  fluo¬ 
rescence  of  formaldehyde  in  an  optical  motored  en¬ 
gine  [15],  The  adiabatic  core  gas  model  [16]  is  based 
on  the  assumption  of  an  adiabatic  core  zone,  sur¬ 
rounded  by  a  boundary  zone  that  cools  by  heat  trans¬ 
fer.  Pressure  is  assumed  homogeneous  in  the  cham¬ 
ber,  but  not  the  temperature.  The  core  gas  volume 
shrinks  with  time  and  vanishes  when  the  reaction 
time  becomes  sufficiently  long.  Temperature  at 
TDC  was  measured  in  the  center  of  the  rapid  com¬ 
pression  machine  using  laser  Rayleigh  scattering 
[17]  and  agrees  well  with  theoretical  adiabatic  core 
gas  temperatures  [18]  calculated  from  the  initial 
temperature,  pressure  at  TDC,  and  temperature-de- 
pendent  specific  heats  of  the  reactants.  Temperature 
measurements  by  thermocouple  show  that  the  core 
gas  occupied  about  80%  of  the  chamber  [19],  Re¬ 
duced  temperature  inhomogeneities  have  been  ob¬ 
served  before  ignition  and  can  be  explained  by  the 
negative  temperature  coefficient  of  reaction  through 
its  feedback  action  on  the  heat  release. 

Delay  times  for  the  first-  (cool  flame)  and  second- 
stage  ignition  of  n-pentane,  isopentane,  and  neopen¬ 
tane  were  deduced  from  pressure  and  light  emission 
traces.  Reactive  mixtures  consisted  of  fuel/oxygen/ 
diluent  in  the  ratio  1/8/30.08;  when  the  diluent  was 
N2,  the  mixture  was  stoichiometric  pentane/air.  To 
reach  compression  temperatures  over  the  range 
640-900  K,  part  or  all  of  the  nitrogen  was  replaced 
by  argon  or  carbon  dioxide.  For  each  fuel,  12  stoi¬ 
chiometric  mixtures  were  studied.  At  least  three 
identical  experiments  were  conducted  on  each  mix¬ 
ture  with  initial  pressures  of  300  and  400  torr.  Ex¬ 
perimental  ignition  delay  times  are  plotted  versus 
TDC  core  gas  temperature  in  Figs.  1-3.  The  autoig¬ 
nition  behavior  can  be  characterized  by  the  following 
quantities:  the  lowest  temperature  at  which  autoig¬ 
nition  takes  place  (Tj),  the  extent  of  the  negative 
temperature  coefficient  (NTC)  zone  (Tmin,  Tmax),  the 
minimum  and  maximum  ignition  delay  times  within 


Fig.  1.  Delay  times  of  cool  flame  and  total  ignition  for 
stoichiometric  mixtures  of  isopentane/oxygen/inert  at  ini¬ 
tial  pressures  of  400  (a)  and  300  (b)  torr.  Open  (cool  flame) 
and  filled  (total  ignition)  symbols,  experiment;  T,  core  gas 
temperature  at  TDC.  Solid  line,  simulation  including  the 
compression  phase;  dashed  line,  simulation  without  the 
compression  phase;  T,  core  gas  temperature  at  TDC. 

that  zone  (fmin,  fmax),  and  the  extent  of  the  cool  flame 
(CF)  zone  (TCFmin,  TCFmax).  These  values  are  re¬ 
ported  in  Table  1. 

Isopentane  (research  octane  number  [RON]  = 
92.3)  is  the  least-reactive  pentane  isomer,  and  its  ig¬ 
nition  delay  times  in  the  NTC  zone  are  the  longest 
(Table  1).  Ignition  was  not  observed  below  690  K;  a 
clear  negative  temperature  coefficient  of  the  reac¬ 
tion  rate  was  seen  between  730  and  820  K,  and  the 
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Fig.  2.  Delay  times  of  cool  flame  and  total  ignition  for  Fig.  3.  Delay  times  of  cool  flame  and  total  ignition  for 
stoichiometric  mixtures  of  neopentane/oxygen/inert  at  ini-  stoichiometric  mixtures  of  n  -pentane/oxygen/inert  at  initial 
tial  pressures  of  400  (a)  and  300  (b)  torr.  Open  (cool  flame)  pressures  of  400  (a)  and  300  (b)  torr.  Open  (cool  flame) 
and  filled  (total  ignition)  symbols,  experiment;  T,  core  gas  and  filled  (total  ignition)  symbols,  experiment;  T,  core  gas 
temperature  at  TDC.  Solid  line,  simulation  including  the  temperature  at  TDC.  Solid  line,  simulation  including  the 
compression  phase;  dashed  line,  simulation  without  the  compression  phase;  dashed  line,  simulation  without  the 
compression  phase;  T,  core  gas  temperature  at  TDC.  compression  phase;  T,  core  gas  temperature  at  TDC. 


temperature  range  of  700-755  K  in  which  cool 
flames  were  observed  is  narrow. 

n-Pentane  (Fig.  3)  is  the  most  reactive  isomer.  Its 
low  octane  number  (RON  =  61.7)  is  consistent  with 
a  low  ignition  limit  7\  of  650  K  and  the  shortest 
ignition  delay  times  in  the  NTC  zone.  The  NTC  and 
the  cool  flame  zone  extend  over  wide  ranges  of  tem¬ 
perature:  755-845  K  and  665-845  K,  respectively,  at 
400  torr. 


The  ignition  limit  for  neopentane  is  about  670  K, 
intermediate  between  n-pentane  and  isopentane, 
consistent  with  its  octane  number  (RON  =  85.5). 
However,  some  peculiarities  should  be  noted.  The 
ignition  delay  decreases  slowly  with  temperature  af¬ 
ter  about  770  K  at  both  initial  pressures  (Fig.  2),  and 
a  negative  temperature  dependence  is  not  observed. 
Experiments  at  200  torr  (not  presented  here) 
showed  a  weak  NTC  behavior  between  770  and  850 
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TABLE  1 

Autoignition  criteria  for  stoichiometric  hydrocarbon/“air”  mixtures 


Hydrocarbon 

Isopentane 

Neopentane 

n-Pentane 

RON 

92.3 

85.5 

61.7 

p0/torr 

300 

400 

300  400 

300 

400 

IV  K 

700 

690 

670  665 

650 

650 

NTC  Zone 

Slowing  Down  Zone 

NTC  Zone 

Tmln/K 

730 

730 

770  775 

755 

755 

T  /K 

1  max'  lv 

820 

820 

880  890 

850 

845 

tmm/ms 

75 

34 

21.5  5 

13 

6 

fmax/mS 

109 

46 

23.5  11 

18 

8 

CF  Zone 

CF  Zone 

CF  Zone 

TcFmin/K 

700 

700 

685  665 

670 

665 

^CFmax/K 

755 

755 

880  890 

840 

845 

Fig.  4.  Pressure  history  of  a  non-reactive  mixture  n-pen- 
tane/N2/Ar  =  0.0256/0.2048/0.7696.  Symbols  are  exper¬ 
imental  values;  curve  is  the  simulation;  initial  conditions, 
p0  =  400  torr,  T0  =  355  K. 

K.  In  a  high-pressure  flow-reactor  study  of  neopen¬ 
tane  oxidation,  Wang  et  al.  [20]  observed  NTC  be¬ 
havior  but  noted  it  was  much  less  pronounced  than 
for  n-pentane  under  comparable  conditions. 

Another  interesting  feature  of  neopentane  autoig¬ 
nition  is  the  persistence  of  cool  flame  activity  over 
the  entire  range  of  investigated  temperatures.  More¬ 
over,  for  compression  temperatures  above  770  K,  the 
cool  flame  occurs  during  the  compression  stroke. 
The  measured  pressure  at  TDC  is  therefore  influ¬ 
enced  by  heat  release  from  the  cool  flame  and  results 
in  core  gas  temperatures  about  20  K  greater  for 
these  cases  for  neopentane  than  for  the  other  iso¬ 
mers  without  cool  flame  behavior  during  compres¬ 
sion. 


Modeling  of  Ignition  Delay  Times 

The  chemical  kinetic  reaction  mechanisms  used 
for  these  simulations  are  based  on  a  long  history  of 
mechanism  development.  There  have  been  recent 
modeling  studies  of  n-pentane  [4,12,20],  neopen¬ 
tane  (2,2-dimethyl  propane)  [20-22],  and  isopentane 
(2-methyl  butane)  [21],  These  mechanisms  are  re¬ 
fined  continually,  but  the  three  mechanisms  used 
here  are  at  the  same  level  and  detail  as  that  reported 
[21]  for  neopentane.  The  three  mechanisms  were 
combined  into  a  single,  larger  mechanism  which  is 
too  large  to  include  here  but  is  available  from  the 
authors.  No  mechanism  modifications  were  made  to 
alter  the  computed  results  or  their  level  of  agree¬ 
ment  with  experimental  data. 

The  numerical  model  computed  reactant  evolu¬ 
tion  through  the  compression  stroke  and  the  sub¬ 
sequent  constant-volume  period.  During  compres¬ 
sion,  there  is  considerable  turbulence  and  mixing  of 
the  reactants  in  the  combustion  chamber,  so  heat 
losses  to  combustion  chamber  walls  can  be  pre¬ 
sumed  to  be  distributed  rather  uniformly  throughout 
the  mixture.  Once  the  piston  stopped,  heat  losses 
were  assumed  to  come  from  gases  near  the  chamber 
wall,  leaving  an  adiabatic  core  that  retained  heat  for 
a  longer  period  of  time.  Therefore  heat  losses  were 
included  for  the  compression  stroke  and  neglected 
during  the  constant  volume  period.  This  assumption 
should  be  reasonably  valid  except  for  long  autoigni¬ 
tion  times  after  the  end  of  compression,  when  heat 
transfer  from  the  core  renders  it  no  longer  adiabatic. 
To  include  the  effects  of  heat  losses  to  the  chamber 
walls,  a  non-reactive  experiment  was  carried  out,  re¬ 
placing  oxygen  with  additional  nitrogen.  A  time-de- 
pendent,  distributed  heat-transfer  coefficient  was 
then  determined  to  be  used  during  the  compression 
stroke  so  the  computed  pressure  history  matched 
the  experimental  results,  as  seen  in  Fig.  4. 
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Fig.  5.  Temperature  history  for 
pentane  isomers  at  the  same  com¬ 
pression  temperature. 


In  another  series  of  model  computations,  the  pres¬ 
sure  measured  at  TDC  and  the  resulting  core  gas 
temperature  were  used  as  initial  conditions,  without 
simulating  the  compression  stroke  and  considering 
the  system  as  adiabatic.  The  results  shown  as  dashed 
lines  in  Figs.  1-3  are  not  veiy  different  from  those 
in  which  simulations  of  reactions  during  the  com¬ 
pression  stroke  were  included.  Therefore  the  con¬ 
cept  of  a  core  gas  at  TDC  is  a  satisfactoiy  approxi¬ 
mation  for  modeling  autoignition  within  the 
operating  conditions  of  the  Lille  RCM,  except  in 
those  cases  where  long  ignition  delay  times  affect  the 
core  temperature. 

Results 

Computed  results  are  shown  as  curves  in  Figs.  1— 
3,  again  showing  the  time  of  onset  of  the  first-stage 
ignition  and  the  time  of  the  final  ignition.  Results 
computed  when  the  compression  period  was  not 
modeled  are  shown  as  dashed  curves. 

Overall  agreement  between  computed  and  exper¬ 
imental  results  is  generally  veiy  good  for  both  total 
ignition  delay  and  the  time  of  the  first-stage  ignition. 
Only  for  the  lower-pressure  experiments  with  iso¬ 
pentane  are  differences  quantitatively  significant, 
with  computations  indicating  ignition  at  about  60 
ms,  while  experiments  showed  ignition  at  about  100 
ms.  At  these  long  experimental  ignition  delaytimes, 
the  assumption  in  the  model  of  adiabatic  conditions 
may  be  inaccurate.  To  test  this  hypothesis,  computed 
constant-volume  portions  of  these  histories  were  re¬ 
peated,  introducing  heat  loss  with  the  same  rate  as 
used  during  the  compression  stroke.  With  the  same 
heat-loss  rate,  these  simulations  failed  to  ignite. 
Even  with  smaller  rates  of  heat  transfer,  these 
mixtures  either  failed  to  ignite,  or  ignited  at  very 
long  elapsed  times.  It  is  apparent  that  differences 
between  computed  and  experimental  ignition  delay 
times  for  isopentane  are  due  to  uncertainties  in  the 
heat-transfer  treatment  in  the  model  and  that  the 


adiabatic  core  in  the  experimental  problem  is  af¬ 
fected  by  heat  transfer  for  reaction  times  of  more 
than  about  40  ms.  In  all  cases  for  the  other  isomers, 
ignition  was  rapid  and  vigorous  in  both  the  model 
and  experiment,  and  the  adiabatic  core  was  insen¬ 
sitive  to  heat  losses. 

The  relative  behavior  of  the  isomers  of  pentane  is 
illustrated  in  Fig.  5,  showing  computed  tempera¬ 
tures  for  each  isomer  at  a  compression  temperature 
of  757  K  and  initial  pressure  of  400  torr.  The  com¬ 
pression  histories  of  the  three  cases  are  virtually 
identical,  but  the  constant-volume  ignitions  are  quite 
different.  All  three  show  a  distinct  first-stage  igni¬ 
tion,  but  the  time  delays  for  that  first  ignition  are 
very  short  (2.4  ms)  for  n-pentane,  longer  for  neo¬ 
pentane  (6  ms),  and  much  longer  for  isopentane  (7.5 
ms).  Temperature  increase  (AT  =  163  K,  106  K,  70 
K,  respectively)  and  fuel  consumption  (/(fuel  = 
41%,  21%,  and  20%,  respectively)  during  first-stage 
ignition  for  n-pentane,  neopentane,  and  isopentane 
are  quite  different,  as  is  the  duration  of  the  first- 
stage  ignition  (0.7  ms,  3  ms,  2.6  ms,  respectively). 
Onset  of  the  second  stage  occurs  in  the  same  order 
as  expected  from  octane  ratings  for  these  fuels. 

For  each  pentane  isomer,  low-temperature  oxi¬ 
dation  takes  place  via  alkylperoxy  radical  isomeriza¬ 
tion  [e.g.,  23,24],  At  these  temperatures,  the  main 
reactions  for  pentane  isomers  are  H-atom  abstrac¬ 
tion  from  the  fuel,  primarily  by  OH 

CsH12  +  OH  =  CsHn  +  H20  (1) 

H-atom  abstraction  by  OII  varies  weakly  with  the 
type  of  C-H  bond  being  broken,  with  primary  C-H 
bonds  being  somewhat  stronger  than  secondary  C- 
H  bonds,  and  tertiary  bonds  being  weakest,  so  II- 
atom  abstraction  will  produce  each  type  of  alkyl  rad¬ 
ical  possible  for  each  isomer.  Processes  leading  to 
first-stage  ignition  begin  by  addition  of  molecular  ox¬ 
ygen  to  the  pentyl  radicals: 
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C5H11  +  02  —  C5Hn02  (2) 

Pentylperoxy  radicals  then  isomerize  via  internal  ab¬ 
straction  of  H  atoms.  Rates  of  these  isomerizations 
depend  primarily  on  two  major  elements,  the  type 
of  C-H  bond  being  broken  and  the  number  of  atoms 
in  the  transition  state  ring  through  which  the  H  atom 
is  abstracted.  The  activation  energy  Eiso  for  isomer¬ 
ization  can  be  written  [24]  as  the  sum  of  ring  strain 
energy,  Ers;  the  bond  energy  for  the  bond  broken, 
Ebond’  and  the  endothermicity,  AE ,  of  the  overall  re¬ 
action. 

EiSO  =  Ers  +  Ebond  +  dE  (3) 

The  most  reactive  pentane  isomer  is  n-pentane, 

II  H  H  H  H 

HC-C-C-C-CH  (4) 

H  H  H  H  H 

This  molecule  can  be  rewritten  to  indicate  different, 
logically  distinct  H  atoms: 

1  2  3  2  1 

1C-C-C-C-C1  (5) 

1  2  3  2  1 

Isomerization  of  pentylperoxy  radicals  in  this  fuel  is 
facilitated  by  the  ready  availability  of  secondary  site 
H  atoms  at  sites  (sites  2  and  3  in  the  model  above) 
that  involve  six-membered  transition  state  rings,  be¬ 
cause  Ers  is  a  minimum  for  a  six-membered  transi¬ 
tion  state  ring,  while  Ebond  Is  smaller  for  breaking  a 
secondary  bond  than  a  primary  bond.  Model  com¬ 
putations  show  that  the  major  products  of  1-pentyl- 
peroxy  and  2-pentylperoxy  radical  isomerizations  are 

H 

O 

OH*  1 1  H 

IIC-C-C-C-CH  (6) 

H  H  H  H  H 

and 

H 

O 

II  O  H  •  II 

IIC-C-C-C-CH  (7) 

HHHHH 

both  favored  because  they  involve  six-membered 
transition-state  rings  and  weakly  bound  secondary 
site  H  atoms,  that  is,  (1,5s)  reactions.  The  rate  of  the 
comparable  isomerization  of  3-pentylperoxy  radicals 
is  relatively  inhibited  because,  although  a  six-mem¬ 
bered  transition-state  ring  structure  is  attainable,  it 
involves  internal  abstraction  of  an  H  atom  from  a 
primary  (l,5p)  process,  rather  than  a  secondary  site. 
Subsequent  reaction  steps  leading  to  chain  branch¬ 
ing  and  rapid  reaction  involve  addition  of  another  02 


molecule  to  the  QOOH  radical  site,  followed  by  a 
series  of  further  internal  H-atom  abstractions,  fol¬ 
lowing  the  principles  for  activation  energy  barriers 
noted  above.  These  final  species  eventually  decom¬ 
pose  to  produce  three  or  more  radical  species,  pro¬ 
viding  chain  branching. 

Neopentane  is  unique  in  that  all  of  its  II  atoms 
involve  primary  C-H  bonds: 

H 

HCH 
H  I  H 

HC-C-CH  (8) 

H  I  H 

HCH 

H 

Following  H-atom  abstraction  and  addition  of  mo¬ 
lecular  oxygen,  there  are  nine  H  atoms  available  for 
internal  abstraction,  each  of  which  is  relatively 
strongly  bound  at  primary  sites,  but  all  have  the  ad¬ 
vantage  of  low- activation  barrier  six-membered  tran¬ 
sition-state  ring  locations.  That  is,  the  pentylperoxy 
radical  has  the  structure 

H 

HCH 

H  !  H 

HC-C-COO  (9) 

H  !  H 
HCH 

H 

Although  the  energy  barrier  to  this  (l,5p)  H-atom 
abstraction  is  1.75  kcal/mol  higher  than  the  (1,5s) 
transfers  in  n-pentane,  there  are  nine  H  atoms  avail¬ 
able,  which  facilitates  this  reaction  step  and  allows 
the  reaction  sequence  to  proceed  relatively  rapidly. 
The  absence  of  an  NTC  regime  for  neopentane  is 
due  in  large  measure  to  the  fact  that  there  is  no 
conjugate  alkene  that  can  be  formed  by  neopentane 
because  of  its  highly  symmetric  and  completely  pri¬ 
mary  C-H  structure. 

Isopentane  has  the  structure 

H 

HCH 

H  !  H  H 

HC-C-C-CH  (10) 

H  H  H  H 

The  most  easily  abstracted  H  atom  is  at  the  tertiary 
site.  However,  subsequent  internal  H-atom  abstrac¬ 
tion  is  inhibited  relative  to  neopentane  and  n-pen¬ 
tane.  In  neopentane,  there  are  nine  H  atoms  at  pri¬ 
mary  sites  with  six-membered  transition-state  rings, 
while  in  isopentane,  there  are  only  three  (l,5p)  H 
atoms  available.  All  other  internal  H-atom  abstrac¬ 
tion  reactions  in  isopentane  are  more  inhibited  be¬ 
cause  they  involve  five-membered  transition-state 
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rings,  either  (l,4p)  or  (1,4s)  reactions,  or  the  initial 
H-atom  abstraction  is  inhibited  because  it  involves 
the  more  strongly  bound  primary  C-H  bonds. 
Therefore,  the  overall  rate  of  chain  branching  in  iso¬ 
pentane  is  lower  than  that  in  neopentane  by  at  least 
a  factor  of  3,  and  its  ignition  rate  is  considerably 
slower  than  that  of  neopentane. 

The  first-stage  ignition  ends  when  the  tempera¬ 
ture  has  increased  enough  to  reverse  the  addition 
reactions  of  molecular  oxygen  to  alkyl  and  hydro- 
peroxyalkyl  radicals.  The  alkylperoxy  radicals  then 
decompose  via  lower  activation  energy  chain  prop¬ 
agation  pathways  to  produce  olefin  +  OH  and  epox¬ 
ide  +  H02  species,  effectively  shutting  off  the 
chain-branching  reaction  pathways,  and  reaction  re¬ 
mains  slow  until  higher  temperatures  activate  new, 
alternative  reaction  pathways. 

For  till  three  isomers  of  pentane,  the  second  stage 
or  hot  ignition  has  exactly  the  same  kinetic  source. 
As  seen  in  Fig.  5,  all  three  fuels  ignite  at  approxi¬ 
mately  the  same  temperature  of  about  950  K.  For 
the  pressures  of  this  study,  this  is  the  temperature  at 
which  hydrogen  peroxide  (H202)  decomposes  at  a 
significant  rate.  Low-temperature  reactions  produce 
an  increasing  amount  of  H202  that  is  stable  until  its 
decomposition  temperature  is  reached;  since  this  re¬ 
action  is  part  of  the  H2/02  reaction  submechanism, 
it  is  independent  of  the  specific  hydrocarbon  fuel 
being  studied.  Past  kinetic  modeling  studies  have 
demonstrated  that  this  same  decomposition  reaction 
is  responsible  for  the  onset  of  knock  in  spark-ignition 
engines  [21,25,26],  ignition  in  diesel  engines  [26,27], 
and  ignition  under  homogeneous  charge  compres¬ 
sion  ignition  (HCCI)  conditions  [26,28]. 

Kinetic  modeling  shows  that  the  isomers  of  pen¬ 
tane  ignite  at  different  times  primarily  because  their 
different  first-stage  ignitions  cause  them  to  arrive  at 
the  hot  ignition  temperature  at  different  times.  The 
most  reactive  isomer,  n-pentane,  has  the  most  vig¬ 
orous  first-stage  ignition  and  the  greatest  associated 
amount  of  temperature  increase,  while  neopentane 
has  a  slower  first-stage  ignition  with  less  temperature 
increase.  Isopentane  has  the  slowest  first-stage  ig¬ 
nition  and  the  least  heat  release.  These  differences 
in  first-stage  ignition  properties  are  attributable  to 
differences  in  their  molecular  structure  and  the 
types  of  C-H  bonds  in  these  molecules.  Following 
the  end  of  the  first-stage  ignition,  n-pentane  is  clos¬ 
est  to  the  onset  of  hot  ignition  determined  by  H202 
decomposition. 

It  is  clear  that  low- temperature,  alkylperoxy  radi¬ 
cal  isomerization  plays  a  major  role  in  many  practi¬ 
cal,  applied  ignition  situations,  including  engine 
knock  and  diesel  ignition.  In  addition,  the  overall 
success  of  the  kinetic  model  to  describe  the  experi¬ 
mental  results  indicates  that  the  most  important  ki¬ 
netic  features  of  these  problems  are  captured  by 
these  reaction  mechanisms,  in  particular  the  influ¬ 


ences  of  fuel  molecular  size  and  structure.  Finally, 
the  value  of  the  RCM  to  examine  these  issues  has 
again  been  confirmed;  in  addition,  this  work  has 
been  able  to  examine  the  sensitivity  of  some  RCM 
experiments  to  details  in  the  heat-transfer  model  and 
spatial  inhomogeneities  in  temperature  in  the  com¬ 
bustion  chamber. 
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COMMENTS 


Eliseo  Ranzi,  Politecnico  di  Milano,  Italy.  Your  low-tem¬ 
perature  oxidation  mechanism  is  very  accurate  in  predict¬ 
ing  the  ignition  delay  times  of  the  different  isomers.  In  my 
experience  with  n-heptane  oxidation  in  RCM,  the  kinetic 
model  overpredicts  the  fuel  conversion  during  the  cool 
flame.  Can  you  explain  how  far  your  model  is  also  able  to 
predict  this  intermediate  conversion? 

Au  thor’s  Reply.  The  experimental  fuel  conversions  of  the 
three  pentane  isomers  during  the  cool  flame  were  not  mea¬ 
sured  in  the  rapid  compression  machine  in  this  work.  The 
ability  of  the  low-temperature  oxidation  mechanism  to  pre¬ 
dict  the  evolution  of  the  chemical  species  remains  to  be 
tested. 

• 

John  Simmie,  National  University  of  Ireland  Galway, 
Ireland.  The  observed  postcompression  decrease  in  pres¬ 
sure  (and  therefore  in  temperature)  in  the  non-reactive 


mixture  in  your  rapid  compression  machine  experiments  is 
incompatible  with  the  neglect  of  heat  losses  in  the  model 
simulation  for  comparable  reactive  mixtures.  Furthermore, 
if  the  simulations  which  do  make  due  allowance  for  heat 
losses  from  the  core  gas  do  not  predict  ignition,  which  does 
occur  in  practice,  then  the  model  is  deficient.  It  is  not  cred¬ 
ible  that  a  model  simulation  without  including  heat  loss 
could  correctly  predict  ignition. 

Author’s  Reply.  The  observed  postcompression  pressure 
decrease  in  the  non-reactive  mixture  experiments  is  not 
related  to  the  core  gas  temperature  but  to  the  mean  gas 
temperature.  The  heat  loss  coefficients  used  in  some  of  the 
simulations  do  not  take  into  account  the  spatial  tempera¬ 
ture  gradient  in  the  combustion  chamber  and  do  not  de¬ 
scribe  the  core  gas  temperatures.  We  agree  with  the  im¬ 
plication  from  Prof.  Simmie  that  the  best  solution  would 
be  a  spatially  variable  computation  that  would  combine  de¬ 
tailed  chemical  kinetics  as  well  as  spatial  heat  transport  to 
the  combustion  chamber  walls. 


